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Abstract: 

In this work, we calculate the decay rates of the OZI-forbidden processes 
■ J/ip(Y) — > 7T7r,/?7r at the order of the leading-twist distribution amplitude. 

The process of J/ip(T) — > tt + tt~ violates isospin conservation and the ampli- 
tude is explicitly proportional to the isospin violation factor m u — rrid, our 
Ch ' numerical results on their decay rates are consistent with the data. The pro- 

£3 , cess J/Tp(Y) — > pn violates the hadronic helicity conservation and should be 

suppressed, as indicated in literature, its decay rate can only be proportional 
to to 2 , at the order of leading twist. Our theoretical evaluation confirms this 
statement that the theoretical evaluation on r(J/-0(T) — > pir) is almost one 
order smaller than the data unless the model parameters take certain extreme 
C"| . values. It may imply that the sizable branching ratio of J /tp(T) — > pir should 

Oh' be explained by either higher twist contributions or other mechanisms. 
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I. INTRODUCTION 

It is generally believed that the narrowness of the ground states of heavy quarkonia J/ ip and T is due to the so called 
OZI suppression This OZI rule demands that if there are no quark lines connecting the initial and final hadron 
states, the processes are suppressed. At beginning, it seemed to be a phcnomenological principle, however, further 
studies indicate that the suppression may originate from the loop suppression which can be precisely evaluated in 
the framework of perturbative QCD. More than 20 years ago, the OZI-suppressed radiative decays of orthoquarkonia 
was investigated by Korner et al. in perturbative QCD[2|, where reasonable approximations were adopted. Since 
then, technique for calculating loop diagrams has been greatly improved and knowledge on the wavefunctions of 
light mesons is much enriched. Meanwhile more data have been accumulated and the corresponding experimental 
measurements become more preciseJ3,0], all the experimental progress indeed provides us with a possibility to test our 
theoretical framework where the perturbative and non-perturbative effects are factorized and a convolution integral 
over them results in the physical transition amplitude. Following their work, we have also rc-calculatcd the rates 
of J/tp(T) — > 7 + 7T° (77, rf) which are respectively isospin-violated, flavor-SU(3)-violated and fiavor-SU(3)-favored 
processes without any approximations at one-loop level[5[. 

In fact, there may exist other possible mechanisms which also contribute to the concerned processes of J/ip(T) — ► PP 
and VP where P and V stand for pseudoscalar and vector mesons respectively [1, 0, [1], therefore to fully understand 
such reactions, a complete calculation on the OZI-suppressed non-leptonic decay processes is obviously necessary and 
should be possible with our present knowledge. Comparing with the radiative decays, theoretical evaluation of the rate 
of the non-leptonic decays is much more complicated. In the radiative decays, a photon is emitted as a free particle 
escaping away from the reaction and it does not participate in strong interaction. For the non-leptonic decay, the 
two (at least) daughter hadrons tangle together by exchanging gluons, therefore one not only needs to carry out the 
complicated Feynman integrations of four-point and five-point loop functions (i.e. D- and E- functions), but also there 
are more Feynman diagrams than the radiative decays. In this work we obtain the transition amplitude by carefully 
calculating the loop integrations. Following the standard procedure^, one can reduce the 5-point loop functions into 
4-point and 3-point loop functions which arc then evaluated in terms of the program "LoopTools" [hJ E3 . Moreover, 
one needs to carefully handle the color factors whereas they are much simpler in the radiative decays. 

In this work, we are going to make a full calculation on the OZI-suppressed processes of J/i{j(T) — > irir and 
J/ip(¥) — ► pir at the order of leading twist. 

The reason to only consider n^' and p*' as the produced pseudoscalar and vector mesons is following. The 
processes are non-leptonic decays, at least three hadrons are involved and to theoretically evaluate the rates, one not 
only needs to calculate the complicated loop integrations at quark-gluon level, but also have to deal with the hadronic 
matrix elements which are fully governed by the non-perturbative QCD effects. However, at present, a completely 
reliable way to calculate the non-perturbative QCD effects is lacking, so that some phenomenological models must be 
invoked. In the decays of J/ip(T) — > tttt, pir, the product mesons arc light and can be nicely described in terms of the 
light-cone distribution amplitudes. Since ir and p are composed of only u, d and u, d whose masses are approximately 
equal, due to the obvious symmetry, the distribution functions are more symmetric and reliable, at least for the leading 
twist order. By contraries, for the distribution functions of K(K*), r\ and rf, the produced mesons are composed of 
constituents u(d) and s quarks which have a large difference in mass, thus one would expect larger uncertainties in 
the evaluation of hadronic matrix elements. Therefore, in this work, these final states are not concerned. 

A simple analysis indicates that J/ip(T) — > irir is an isospin-violating process. Namely the pions are treated as 
identical particles once the isospin symmetry is adopted in the analysis. Concretely, by the conservation of angular 
momentum, the two pions are in the p-wave state, since pions are identical bosons, the wave function of the two-pion 
system must be totally symmetric, so that the isospin of the system should be 1 as 

i)|i, -1) - |i, 1)) = -^(k + }IO - IOk + »- 

That requires that the process of J/ip — > 7r°7r° is strictly forbidden. The isospin violation effects are expressed in the 
mass difference of u and d quarks which appears at the loop calculations, and the factor m u — will be explicitly 
shown in the expressions derived at the quark level. Even though, we only consider the leading twist contribution of 
the light-cone wave functions which arc independent of the quark masses, we still count in the mass splitting which 
results in the isospin violation. Turn to the processes J/ip(T) — > pir, in contrast with the irir case, p and ir are not 
identical particles, therefore the anti-symmetry requirement which enforces the two-pion system to be in isospin 1 
state, is dismissed, so that the pir system can be in isospin state and it guarantees the iso-spin conservation for 
the decay process J/ijj(T) — > p Q ir°. This observation seems to demand that the branching ratio of J/tp(T) — > pir 
should be much larger than the isospin violating process J/ip(T) — > irir and the data indeed support this statement. 
Moreover, in J/ip(T) — > pir, the isospin state of pir is dominant, so that the branching ratios of J/ip(Y) — > p°ir° 
and J/ip(Y) — > p + ir~ + p~ir + roughly retain a relation of 1:2. However, as indicated in Rcfs.fll, [Hj], such processes 
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violate the hadronic helicity conservation because gluons and photon do not carry hadronic helicitics. A non-zero 
theoretical prediction on the rate at the order of leading twist must come from a violation of the hadronic helicity 
conservation. It is indicated that such a violation is proportional to the light quark mass, therefore one can expect that 
the directly calculated OZI suppressed amplitude should be proportional to rri^(q = u,d). Our calculation confirms 
this mechanism (see the text for details). 

The data[lil tell us an opposite conclusion that the helicity-violated process J/V>(T) — ► pit has sizable branching 
ratio and almost is one of the dominant modes in J/ip and T decays. The discrepancy should be explained, some 
suggestions that the next-to-lcading twist contribution, higher Fock states and other mechanisms such as the hadronic 
loop and glueball intermediate states etc. are taken into account, are proposed. 

In this work, we only concern the OZI-supprcssed processes for J/ip(T) — > tttt, pir, and will explicitly demonstrate 
the isospin conserving and violating effects and the helicity violation effects in the formulation. 

As indicated above, to evaluate the hadronic matrix elements, one has to deal with a convolution integrals over 
the distribution amplitudes of the concerned hadrons. Because J/ip and T contain two heavy constituents, their 
bound-state effects can be simply expressed in terms of the wave functions at origin which can be easily obtained 
from the data of their leptonic decays. The distribution functions of the two produced light hadrons might cause 
uncertainties, even though as indicated above, for n and p mesons, they can be reduced to minimum. It is interesting 
to note that the OZI-suppressed process for J/tp(Y) — > pir is forbidden by the hadronic helicity conservation at the 
leading twist, if the quark mass is neglected at the loop calculations. However, it is not zero and the transition 
amplitudes of such processes must be proportional to m^. In this work we only consider the contribution from the 
leading twist distribution amplitudes of the mesons and show that as m q — ► 0, the amplitudes would approaches zero, 
in other words, we confirm the statement that the hadronic helicity conservation forbids the process J/ip(T) — > pir at 
the leading twist if quark mass is neglected. 

Following the literature, we can trust the calculations to a relatively accurate level. A rough numerical estimate 
by changing the input parameters in the distribution functions and its forms given in literature, shows that the error 
can be of order of a few tens percents. 

After this introduction, we give all the formulas where we carry out the four- and five-point Feynamn integrals to 
obtain the hard-scattering amplitude at the quark level. The isospin violation factor m u — explicitly shows up 
for J/ip(T) — ► 7T7T, and for the helicity-violated J/"4>{T) — > p7r the amplitude is also proportional to the light-quark 
masses, then one needs to convolute the hard kernel with the initial and final states, and the convolution integration 
results in the physical transition amplitude in Sec. II. In Sec. Ill, we carefully analyze the infrared behavior in 
J/^(T) — ► 7T7T, pir and convince ourselves that all Feynman diagrams are infrared-safe when the end-point behaviors 
of the wave functions are considered. In Sec. IV, we make a numerical evaluation of the decay rates of J/ifj(T) — > ir + ir~ 
and p ± 7T zfl and some necessary input parameters are explicitly given. The last section is devoted to a simple discussion 
on the uncertainties in our calculation and possible contributions to these processes from other mechanisms and then 
draw our conclusion. Some tedious details are collected in the appendices. 



II. THEORETICAL CALCULATION ON THE RATES OF J/ip(T) -> tttt^tt 

In this work, without invoking the so-called weak-binding approximation which was adopted in literature Q, we 
explicitly keep the masses of the heavy and light quarks at the concerned propagators, when derive the transition 
amplitudes. The amplitude is written as 

A = H ® $jy^(<r) ® $Pi ® ®Pz 

H = C®H (1) 

where the factors C, H, ^j/^(T),p lt p 2 are the color factor, hard kernel and distribution amplitudes of mesons, re- 
spectively. And the labels P\,P2 denote the two produced mesons in the final state. Indeed, here the perturbative 
and non-pcrturbative parts arc factorized and a convolution integral would associate them to result in the physical 
amplitude. The detailed expressions of the hard kernels are given in Appendix A. 

Below, in Fig. 1, we present the relevant Feynman diagrams. In these figures, we only explicitly draw the typical 
diagrams. There exit their topologically deformed diagrams which are obtained by exchanging the connections of the 
gluon-lincs in the loop to the light-quark-gluon vertices, namely the two gluon-lines cross with each other. We do 
not explicitly show them in Fig. 1 just for simplicity. But definitely, in our derivation the contributions from those 
diagrams are included. 

The amplitude of J/i/)(T) — > P\?2 can be divided into three categories which correspond to Fig. [1] (lal) and (la2), 
(lbl) and (lb2), (lcl) and (lc2) respectively. To make the text succinct we collect the detailed expressions of the 
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FIG. 1: The Feynman diagrams of J/tp(T) — > P1P2. Our calculation also includes such diagrams which are topologically 
deformed from that shown above by exchanging the connections of the gluon-lines in the loops to the gluon-light-quark vertices, 
where the two gluon-lines cross with each other. 

amplitudes in Appendix B, where the diagrams with gluon-lines in the loops crossing each other are labelled as (2al), 
(2a2), (2bl), (2b2), (2cl) and (2c2) respectively. 

Generally, the transition amplitude for J/ip(Y) — ► 7r + 7r~, p + ir~, p~ir + can be written as: 

i 
i 

where summing over i means including all diagrams listed above and their topologically deformed diagrams which 
were depicted above. For tt + ti~ final states one possible setting is that Pi, P2 correspond to ir + , ir~ respectively, 
and another possibility of interchanging tt + and 7r~ is also included in the sum. m q is either m u or rrid, which 
respectively exist in different settings and their contributions should be summed in the final amplitude. For p + ir~ or 
p~ir + final states interchanging Pi, P 2 would induce obvious differences and therefore we use two new labels "a" and 
"b" to distinguish the two different settings. Here we do not need to calculate the rate of J/ifi(T) — > p°Tr° because as 
discussed above, it is an isospin conserving process and the Clebsch-Gordan coefficients in 1=0 state determines the 
ratio of T(J/i/j(Y) — > p°Tr°)/r(J/t(j(T) — > pn) and it should be close to 1/3, both the data and the analysis according 
to the topology of our diagrams shown in Fig. [T] confirm it, even though we only consider the contributions from the 
leading twist distribution amplitudes of the mesons. 

In the quark picture, hadrons are made of valence quarks whose momenta-distributions are described by appro- 
priate distribution functions. The leading-twist distribution amplitudes of J/ip(T) is usually defined through the 
correlator[l5j: 

^fj/i>rnj/^j/^j/^\\(u) + -a 11 v «/,//</, (ej/,^'?v - Ej/^'TV^J/iM") j 7 (3) 

where £j/^ and fj/^ are the polarization vector and decay constant of J/ip respectively, and u = 1 — u. 0ii and <j)±_ 
are the leading-twist distribution functions corresponding to the longitudinally and transversely polarized mesons, 
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respectively, by the definition given in literature [15(. For the case of T one only needs to replace all the symbols of 
charm c into bottom b (of course as well as the concerned parameters). The leading-twist distribution amplitude of 
7r is usually defined through the correlator [l6j: 

WMM^) = i 5 -^ f du^v'-v^-* {^ 5 J>(u)} 0a , (4) 

where f n is the decay constant of pion. And the leading-twist distribution amplitude of p is usually defined through 
the correlator [15j: 

(p(p')\q\ a (y)qi (x)\O) = h-j\ ue ^'-v+^'- x 
f P m p f p 4> p \\{u) - \^ v 'if p {e^, V ' vl - e* pu ,p',,)4>Mu)) (5) 



0a 



where e p and f p , fj are the polarization vector and decay constant of p respectively. 



III. ANALYSIS ON THE INFRARED BEHAVIORS IN J/ip(T) -> PiP 2 

It is definitely demanded that a reasonable theoretical prediction on any practical process must be infrared safe, 
namely the infrared divergence must be exactly cancels if it exists or properly dealt with at the end of the calculation 
which corresponds to real physical measurable quantities, such as the decay width and cross section. In this work, 
we will explicitly show that in the framework of pcrturbative QCD, the infrared behavior of each individual Feynman 
diagram shown in Fig. [I] is benign, even though at first glimpse it seems to be divergent. 

There are several typical Feynman diagrams shown in Fig. [TJ We take the amplitude of Fig. Q](la2) for J/tp — > P1P2 
as an example to analyze the infrared behavior. 

Its contribution to the transition amplitude reads 

A la2 (m q ,u,v) = H la2 (m q ,u,v)®<S> Pl (u)(Z><f>p 2 (v)., (6) 
and H la2 (m q ,u,v) is given in Appendix A. The concerned factor of the amplitude is 



k 2 (k +pi + p 6 ) 2 [{k +Pi) 2 - m q ][(k+pi - p 3 - p 5 ) 2 - m 2 Q ](p d + p 5 ) 2 [(pi - P3 - Ps) 2 - m 2 Q ] ' 

Firstly, if there exists an infrared divergence in the loop integration, it must come from the kinematic region k — > 
and the end-points of the distribution functions, therefore one only needs to analyze two cases: (1) k — > 0, (2) at 
end-points. 

To show the infrared behavior of the amplitude after integrating out the loop function in the case (I), we may fix 
the external momenta of the quarks and antiquarks in final states by a special choice ps — P4 — \p Pl , P5 = P& = \p P2 
to avoid possible endpoint divergence. Looking at the expression ([7]), as m q 7^ 0, the dangerous part is proportional to 
\/k 2 which is finite after integrating over the loop momentum d 4 k, so that in this case there is no infrared divergence 
coming from the loop integration. 

Secondly, in case (2), when the momentum of one quark(antiquark) in each of the final mesons is close to its 
endpoint, for example P4, pe — ► 0, while the other quark (antiquark) takes almost all the momentum of the meson. 
It is observed that the factor ^jj^rp does not contribute a divergence. The dangerous term comes from the factor 

k 2 (k+p4 +po) 2 at the denominator as k — > and P4, p^ — > 0, which seems to cause a logarithmic divergence. However 
as we convolute the amplitude with the distribution functions of the two mesons whose distribution functions linearly 
approach to 0(</>(w), <f>(v) — > 0) at the end-points, because it turns to zero fasters than the logarithmically divergent 
factor, i.e. linin^o win it — > 0, the infrared behavior is safe. 

Finally, when ^3, p§ — > 0, the loop integration does not produce any divergence by the same sake of the first 
case. It is noted, that the factor ( P3 + P5 )2 0( M ) < /'( t ') is finite, but there exists a subtlety. Namely in general the limit 
depends on the ways how x,y approach to more or less. For J/ijj(T) — > itit there are two possible settings for each 
diagram, namely, an interchange «-> tt^ brings one setting to another, and their contributions have an opposite 
sign due to the SU(2) symmetry and cancels each other (obviously, for the finite term, their contributions cannot 
cancel each other due to the SU(2) breaking i.e. m u ^ m<j). Thus the dependence on the order of limits disappears. 
By contraries, for J/ip(Y) — * pn, there is no such a cancellation, so that even though infrared divergence does not 
exist, the final numerical results somehow depend on the order of taking limits of u and v approaching to zero. The 
strategy we adopt in this work is to set the integration order as we integrate over u and then v and it can be treated 
as a rcgularization scheme similar to that we generally adopt for treating the ultraviolet divergence. 
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IV. NUMERICAL RESULTS 

The input parameters which we are going to use in the numerical computations are [El. Ibil. fl5L UtI [bsj : /,//</> = 551 
MeV, f T = 710 MeV, U = 131 MeV, f p = 198 MeV, /J = 160 MeV, = 3096.87 MeV, m T = 9460.3 

MeV, m n ± = 139.57 MeV, m p ± = 775.5 MeV, a s (m c ) = 0.32, a s {m b ) = 0.21, m c = 1300 MeV, m b = 4500 MeV, 
and the meson distribution functions respectively. For the numerical evaluations, in Eqs.([U [5]), we adopt three 
different distribution functions for the light mesons in the literatures [jj, QJ, HJ EH as <^i(x) = 6a: (1 - x), fa(x) = 
30x 2 (l — x) 2 7 <f)3(x) — ^(1 — 2.t) 2 [1 — (1 — 2x) 2 }, and also let the current quark masses of the u and d types vary 
within a reasonable range. 

Below in Tables I, II, III and IV, we present our numerical results on the decay rates of J/ip — > 7r + 7r~, p + 7r~ + p~n + 
and T — * 7r + 7r _ , p + n~ ± p~ir + respectively. 

TABLE I: Decay widths (T) of J/rp — > ir + -ir~ based on the three distribution functions, 4>\, 02 and 03, respectively 



m„(MeV) 


m d (MeV) 


r(0i)(MeV) 


r(0 2 )(MeV) 


r(0 3 )(MeV) 


exp(MeV) 


2 


4 


4.52 x 10" 5 


2.98 x 10" 5 


2.71 x 10" 4 




3 


4 


1.88 x 10~ 5 


9.35 x 10" 6 


5.67 x 1(T 5 




3 


5 


3.17 x 10~ 5 


2.36 x 10" 5 


1.25 x 10" 4 


(1.37 ±0.21) x 10 -5 


4 


5 


1.03 x 10~ 5 


8.12 x 1(T 6 


4.26 x 10" 5 




4.5 


6 


2.29 x 10" 5 


1.43 x 10" 5 


8.85 x 10" 5 





TABLE II: Decay widths (r) of T — » n + n based on the three distribution functions, 0i, 2 and 03, respectively 



m u (MeV) 


m d (MeV) 


1X00 (MeV) 


r(0 2 )(MeV) 


r(0 3 )(MeV) 


exp(MeV) 


2 


4 


2.79 x 10" 6 


1.24 x 10" 6 


1.13 x 10" 5 




3 


4 


8.16 x 10~ 7 


5.28 x 10" 7 


6.95 x 10" 6 




3 


5 


1.23 x 10~ 6 


9.43 x 10~ 7 


9.6 x 10~ 6 


< 2.7 x 10~ 5 


4 


5 


7.39 x 10~ 7 


2.72 x 10~ 7 


5.11 x 10~ 6 




4.5 


6 


9.78 x 10" 7 


7.5 x 10" 7 


8.93 x 10" 6 





TABLE III: Decay widths (F) of J/tp — > tt + p ± tv p + based on the three distribution functions, 0i, 2 and 03, respectively 



m u (MeV) 


m d (MeV) 


T(0i)(MeV) 


T(0 2 )(MeV) 


r(0 3 )(MeV) 


cxp(MeV) 


2 


2 


1.04 x 10" 4 


7.21 x 10" 5 


5.11 x 10" 4 




3 


3 


2.36 x 10~ 4 


1.6 x 10~ 4 


1.17 x 10~ 3 




4 


4 


4.12 x 10~ 4 


2.9 x 10" 4 


2.08 x 10 -3 


(1.06 ±0.08) x 10" 3 


5 


5 


6.69 x 10~ 4 


4.54 x 10~ 4 


3.38 x 10~ 3 




6 


6 


9.75 x 10" 4 


6.68 x 10 -4 


4.88 x 10" 3 





As discussed above the OZI-suppressed process J/tp(T) — > ir + ir~ is isospin-violated whereas J/0>(T) — * n + p~ ± 
ir~ p + violates the hadronic helicity conservation. It is noted that the theoretically evaluated values on the OZI- 
suppressed processes J/i/j — > tt + tt~ arc slightly larger than the experimental data depending on the parameter choices 
such as m u , md and types of the meson distribution functions, whereas that on J/tjj — > tt + p~ ± tt~ p + are one order 
smaller than the data. It may imply that some other mechanisms may also contribute to the decays and we will 
remark on the results in next section. 
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TABLE IV: Decay widths (r) of T — > tt + p + 7r p + based on the three distribution functions, (j>i, <j>2 and cf>3, respectively 



m„(MeV) 


m d (MeV) 


r(0 1 )(MeV) 


r(^ 2 )(MeV) 


r(0 3 )(MeV) 


cxp(MeV) 


2 


2 


2.23 x 10" 6 


1.56 x 10" 6 


8.54 x 10" 6 




3 


3 


5.04 x 10~ 6 


3.67 x 1(T 6 


1.84 x 10" 5 




4 


4 


8.95 x 10" 6 


6.25 x 10~ 6 


3.4 x 10~ 5 


< 1.08 x 10 -5 


5 


5 


1.42 x 1CT 5 


9.93 x 10~ 6 


5.44 x 10" 5 




6 


6 


1.84 x 10" 5 


1.38 x 10" 5 


7.61 x 10" 5 





V. DISCUSSION AND CONCLUSION 

In this work, we calculate the contributions of the so-called OZI forbidden processes to the decays J/ip(T) — > tttt, pit. 
As we discussed in the introduction, the process J /ipiT) — > irir is an isospin violating reaction, whereas J /ij}(T) — > pir 
is an isospin conserving one, on other aspect, the former one conserves the hadronic helicity whereas the latter one 
violates it. Our numerical results on J/ip(T) — > n + ir~ are reasonably consistent with the data at order of magnitude, 
but the evaluated branching ratio of J/tjj(T) — > pir is obviously smaller than data by one order. 

As shown in Tables I through IV, one can notice that the results deviate from each other in a wider range as one 
adopts different wave functions which all are suggested in literatures, as well as the light-quark masses. 

J/ip(Y) — > 7T7T is an isospin violating process and at the leading twist the OZI-suppressed process which was supposed 
to be the main contribution to the mode of J/ip(T) — * pir violates the hadronic helicity conservation. As well known, 
the source of isospin violation can be either from a photon emission (absorption) and/or quark mass difference, and 
for the helicity violating processes, the decay width is proportional to m 2 ,, so that both of the processes are somehow 
sensitive to the light quark masses and much suppressed. Our formulas explicitly show that as m q — > 0, the decay 
widths for both modes approach zero. This observation confirms the above statements. 

In this work, we only include the contributions from the leading-twist distribution amplitude and our results confirm 
that due to violation of helicity conservation, the theoretical evaluated ratio is one order of magnitude smaller than 
the data. It is also noted from our qualitative analysis that the rate of isospin-violated process J/ip(T) — > irir should 
be proportional to the square of mass difference (m u — irid) 2 , whereas rate of the the hadronic helicity- violating 
process J/ip(T) — > pir is proportional to (m u + m c i) 2 , i.e. it is natural to expect that T(J/ip(T) — > pir) which is 
theoretically estimated in this framework, is a few times larger than T(J/ip(T) — > irir), our numerical results shown 
in Tables I through IV confirm this statement, and if m u — rtid, the estimated T(J/ip(Y) — > irir) is zero, whereas 
r(J/'0(T) — > pir) is not. But this still does not explain the largeness of the branching ratio of J/ijj(T) — ► pir. It is 
indicated in Refs.[ll, HH, the large branching ratio might be due to higher twist contributions. Therefore it seems 
that to correctly evaluate the branching ratio, in principle one needs to include the contributions from higher twist 
distribution amplitudes in the evaluation. 

On other side, besides the contributions from higher twist distribution amplitudes, there may exist other mechanisms 
which may result in larger branching ratios for J/<ip(T) — > pir. As suggested by Suzuki paper [23] and our earlier workQ, 
there can be a contribution from the hadronic loops and by fitting data (in the paper, the contributions from the 
OZI- forbidden processes were not theoretically calculated as we do in this work, but obtained by fitting data), we 
reached two conclusions that if only the two mechanisms contributing, the hadronic loop contribution would have the 
same order of magnitudes as that of the OZI-forbidden processes (definitely including higher twist contributions) and 
secondly the two contributions are destructive. 

Of course, it may not be the end of the story that some authors also suggested a glueball contribution which should 
be added to that from the aforementioned mechanisms @, and then the picture becomes more complicated, because 
we are unable to reliably estimate the glueball mass and phenomenological behaviors so far, unless we can borrow the 
lattice results. Therefore further developments on theory are necessary. 

Uncertainties in our theoretical evaluations come from the input parameters, especially the light quark masses and 
the shapes of the distribution functions while only the leading-twist distribution amplitudes are accounted. One can 
note that the shapes of the wave functions would cause order of magnitude differences. So far, we still cannot really 
rule out any of them, but wait for more accurate data to determine. 

As indicated in the text, we only consider the processes of J/tp(T) — > irir, pir because as there no strange flavor 
gets involved, the wave functions of the produced mesons is simpler and more symmetric. For the processes involving 
such as K(K*), r\, 77', the calculations become more complicated and the results are not much reliable. Therefore we 
postpone our study on such processes in our later works. 

So far, the experimental data are not accurate yet, especially for the measurements on T decays only upper limits 
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are set. However, we are inspired by the promises from the CLEO c and BES III collaborations, as they will provide 
a much larger database on J/ip decays, and more data would be accumulated in the B-meson factories, and then 
we will have concrete numbers about the branching ratios of T — > irn, pir instead of the upper limits set by the 
present experimental measurements. Moreover, the LHCb and future ILC can much enrich our knowledge on hadron 
structure. Conclusion is definite that further work is necessary. 
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Appendix A: The hard-scattering amplitudes H % ' a 13 7 9 al3 ~ ,p {m q ,u,v) 
The hard-scattering amplitude corresponding to Fig. [T](lal) is: 

H lal ' aa 'M'pp'(m q ,u,v) = 

I 7^4^[(-i5.2S7")r— i t~~ j (-ig s Tl^h (-ig s T^ x )} a > a Qi 

J (27t) 4 ft + pi - p 3 - p 5 - rriQ pi - p3 - p5 - tuq 

(Hfe) z p A _ mq (-^i7m)W^ [Hg^7x)]^<g {p3 g {k - ~\ p&)2 



(8) 



The hard-scattering amplitude corresponding to Fig. [I](la2) is: 

H la2 ' aa 'W pp '{m q ,u,v) = 

^Ql^-ig^)—- ^— { - igs Tl^) * {-ig s T^)] a , a Qi 

(2tt) 4 jk+pi- p3 -f>5- m Q pi - p 3 -p5~mQ J 

The hard-scattering amplitude corresponding to Fig. [T](lbl) is: 

H lbl ' aa 'W pp '(m q ,u,v) = 



(9) 



(27r) 4 a ls j)i + P'q - - m Q sr + p'i + fa - fa - m Q 

OS? l(-i9sTMbWwti P ' \H9sT c kw lx) p _ ^ _ - {-i9sT b wmlti )] p . p q? p ^ + ^ )2 g (fc - - — } 

The hard-scattering amplitude corresponding to Fig. [T](lb2) is: 

H lb2 > aa 'f } 0'pp'(m q ,u,v) = 



(10) 



(27r) 4 ft4+P6-p2-rriQ ft + pi + ps - p2 - mQ 

tip [(-igsTtnlu)}^^ [{-i9sT b mwlli ) _^ - ^ — (^g s TUlx)] P > P <ll ( ^ + ^ )2 g (fc _ - pg)2 

The hard-scattering amplitude corresponding to Fig. [T](lcl) is: 

^,[(- Ws7 -)__ * (-^T^)— * (- lgs T r c n x )] a , a Qi 



(11) 



(27TJ 4 ft +pl -pi -p5 - TUQ fi+Pl-TTlQ 

^'[{-WsTnwlu ) — IJ— 7— — — {-igsT^x)]p>f3q2 q[ p \{-igsTi m ^^)] p > p q" l p 

f P3 m q 



-I —I —I 



(P4 + Pb) 2 k 2 (k + pi + p 2 - Pi - PbY 



(12) 
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The hard scattering amplitude corresponding to Fig. [T](lc2) is: 

H lc2 < aa 'M'pP' '(m q ,u,v) = 

lfe r Q J a ,[(-^ 7 ")__ — * (_ iflaT ^y)irrrr 



(27r) 4 ft+pi -p4-p5-mQ jt+pi-mQ J 



t 2p ,[{-ig s n wlx )—- (-%!C„7.)]^/3?2 n /39v[(-^ s ^7M)]p'p?i 



p 

"9 

— 2 — Z —2 



(13) 

(P4 + Ps) 2 fc 2 (fc + Pi + P2 - Pi - Pb) 2 

The hard-scattering amplitude corresponding to the diagram which is topologically deformed from Fig. [T](lal) by 
exchanging the connection of the gluon-lines in the loop to the gluon-light-quark vertices, is 

H 2al ' aa 'W'pp'(m q ,U,v) = 

~ i a >[(-i9.TZT»)- r — 1 t~~ ~7 (-i 9s Tl^)- -p-L —(- igs T^)] a , a Qi 



<%AH9sTLl») f + — Hg.^)W^[Hg.JS,7X)]^ (p3 7 P5 )2 fc? (fc + P 7+Pa) 2 

(14) 

The hard-scattering amplitude corresponding to the topologically deformed diagram from Fig. [T] (la2) is: 

H 2a2 > aa 'W'pp'(m q ,u,v) = 

^A{-i9sW)ir—^ — ~t~ ~i {-ig s T b sr r)n — -n^7 —(-ig s T r c i7 x )} a , a Qi 



(2tt) 4 F + PI - - Z>5 - "IQ Pl-fc-fc- TTlQ 

<?20> [{-WsTi w ^ v ) —g — - 7 (-^T^„7 AI )] /3 ^^ l )3 g I 1 n p ,[(-ig ;s r^ fe 7A)]p' (!) 9ip 



-JS - jfe - w " WJP p "" w LV mfc ,A/JP (p 3 + p 5 ) 2 fc 2 (fc + P4 + Pe) 2 

(15) 

The hard-scattering amplitude corresponding to the deformed diagram from Fig. [1] (lbl) is: 

H 2bl ' aa '^'pp'(m q ,u,v) = 

^Ql[(-ig s T^) \ (-ig s Tl>y») / (- )] Q , Q Q£ 

(27t) 4 p4 + - P2 - mQ jk + pi + pe - p2 - mQ 



(16) 



The hard-scattering amplitude corresponding to the deformed diagram from Fig. [T] (lb2) is: 

H 2b2 ' aa 'M'pp'(m q ,u,v) = 

JpL&ft-ig.TZ'f) , , , ' (-ig.ltV) , , , , / , (-^^- 7 A )]a' Q Q J Q 

(27T) 4 Pi+P6-p2~ m Q fi + Pi+P6-P2- m Q 



jt-p 3 -m q ' w p {Pi + Pe) 2 k 2 (k - p 3 - p 5 )< 

The hard-scattering amplitude corresponding to the deformed diagram from Fig. [T] (lcl) is: 

^Qt,[(- iWs7 -) J— { - iga Tl^)—^ (-ig s T^ x )] ala Qi 

(27r) 4 jk + pi — J>4 — f>5 — niQ jk+pi- itiq J 

<?V [(~^ r n«.7A) p - ^ _ — (-WsT^u^p'nq^q'ip' [(-WsTt m 7^)}p'pqT P 

—i —i —i 



(17) 



(Pi + Pb) 2 fc 2 (fc + Pi + P2 - Pi ~ Pbf 



(18) 
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The hard-scattering amplitude corresponding to the deformed diagram from Fig. [T] (lc2) is: 
H 2c2 < aa 'Wpp'(m q ,u,v) = 



q* ,[(-ig s Tg wlv )——^ — (-^ s T^ l 7A)]/3'/3<Z2V9 J ip'[(-^7 1 l 7 M)]p'p'7i 

—i —i —i 



(P4 +P5) 2 fc 2 (fc +Pl +P2 -P4 - Pa) 2 

Appendix B: The amplitudes A l (m q ,u,v) 
1. For J/V> ->■ PP 

For amplitudes *4 lal and „4 la2 , we have 



(19) 



A iai = c lal H lal (m q ,u,v)^j / ^p 1 ^ P2 

A ia2 = C la2 H la2 (m q7 u,v)<S>j / ^ P1 <S> P2 (20) 



with 



^jlal ^j.^^CLrjnbrj^C^rj^^rj^arj^brj~iC^ 

^jla2 ^^^^a^brj^C^rY^^^brparj^C~j 

ITT 2 1 

H la \m q ,u,v) = -H la2 (m q7 u,v) = -—jg e s ( — f 



'P2 

'pi 



(27T) 43s v 4iV c ' (pa + P5 ) 2 [(Pi -Pa- P 5 ) 2 - m 2 Q ] 
{D (m q ,u,v)[-32mj/^£j/^ ■ p pl p[ ■ p P2 m 2 Q + 96mj/^£ J/4 , ■ p P2 p[ ■ p P1 m 2 Q 
-32mj/^e.j/^ ■ p[p pi ■ p P2 m 2 Q + 32m, 7/ ^£ J/v , • p pl p[ ■ p P2 p 3 - 32m J/4 ,e.j/^ ■ p P2 p[ ■ p pl p 3 
-QAmj/^Ej/^ ■ p pl p[ ■ p' 3 p' 3 ■ p P2 + 64mj/^£j/^ ■ p' 3 p[ ■ p pl p' 3 ■ p T 
-128mj/^ej/^ ■ p P2 p[ ■ p' 3 p' 3 ■ p P1 - 64m J/ ^£ J/ . i/> • p[p 3 ■ p P2 p' 3 ■ p f 
+32m J/ . l pS J/ . l p ■p[p' 3 p pl -p P2 ] 

+D li (m q ,u,v)[-32mj/^e J /^ ■ p P1 p^ P2 m 2 Q + 96m J/4 ,£ J/tl , ■ p^p^rriQ 
-32mj/^^p pl -Pp^rriQ - 32m J/i ,£j/^ ■ p^p^jp'x ■ P 3 - 32m, 7/ ^£ J/v , ■p P2 p^ 1 p / 1 ■ p 3 
+32mj/^ej/^ ■ p P1 p' 3 p[ ■ p P2 + 32mj/^£j/^ ■ p^p^p^ ■ p P2 - 32mj/^£j/^ ■ p P2 p' 3 p[ ■ p P 
+Z2mj/ i ,£j/ ij ■ p 3 p^ 2 pi ■ p P1 + 32m J/i ,£ J/4 , ■ p P1 p^ 2 p 3 - 32mj/^£j/ xjl ■ P P2 P P1 P 3 
-2>2mji^£ji^ ■ p p1 pTp' 3 ■ Pp2 ~ 64m,7/^£j/^, • p P1 p' 3 p' 3 ■ p P2 - 32m J/ ^£ J/ . 4 , ■ p'xP^p^ ■ p P 
+6Am J/i ,£j /i , ■ p' 3 p^ pi p' 3 ■ p P2 + 32m J/i ,£ p J/i> p' 1 ■ p pl p' 3 ■ p P2 - 96m J/i} £ J/i , • p P2 p'fp' 3 ■ p P1 
-128mj /v ,£,7/^ • p P2 p' 3 p 3 ■ p P1 - 32m J / i ,£j/^ ■ p[p^ 2 p 3 ■ p P1 - 32m J / i ,£^ / ^p' 1 ■ p P2 p 3 ■ p f 
-64mj/^ /ifi p 3 ■ p P2 p 3 ■ p pl - 32m J/l> £j/ 1> • p' 3 p?p P2 ■ p P1 + 32mj/^£j/^ ■ p[p' 3 p P2 ■ p P 
+32m J /^ J/ ^p' l -p 3 p P2 -p P1 +32m J/i ,£> 1 J/i> p 3 2 p P2 ■ p pl ] 
+D^(m q ,u,v)[-6Am J / iJ £j/ 4 , ■p P2 p' 3 P pl +64m J/4 ,£ J/4 , ■ p' 3 p^ 2 p pl 
-32m J/4 ,£ J/4 , ■ p pl g^p' 3 ■ p P2 - 96m J/4 ,£ J/4 , ■ p P2 g^p' 3 ■ p P1 - 64m j , N ,p p2 p' 3 ■ p P1 
-32mj/^£j/^ -p' 3 g^p P2 ■ p P1 + 64m.j/^ /4i p 3 p P2 ■ p P1 ]} 

If 1 

D (m q ,u,v) = — I d A k 



DJm q ,u,v) = -r-r f d 4 k 



S-2 I d4k k 2 [{k+pi) 2 _ m 2 q]{k + p4+ p 6 )2 [{k + pi -p 3 ~ p 5 )2 - m 2 Q ] (21) 



k 2 [(k+p A ) 2 


- m 2 ](k - 


h Pi + Pe) 2 [(k + pi -7J3 


-p 5 ) 2 -m 2 Q ] 


k 2 [(k+p 4 ) 2 


— m 2 ](k - 


VPi +P6) 2 [(k +pi - P3 
k^kv 


-P5) 2 - m 2 Q ] 



where p[ = p 4 , p' 3 = pi - p 3 - p 5 . 

For amplitudes _4 lbl and A lb2 , we have 



with 



^161 = c lbl H lbl (m q ,u,v)^j / ^p 1 ^ P2 
A ib2 = c lh2 H lh2 (m q ,u,v)^j / ^p 1 ^ P2 



^<161 r^^^rparpbrpC^r^/rpbrparpC^ 



H (m q ,u,v) = — H (m Q , u, v) = — , o ( Ar ) t^t- ^ ^ 

{D (m q ,u,v)[96mj/^£j/^ -p Px p' x -p P2 mQ - 32m J/i) e J/l p -p P2 p\ ■ p P1 m 2 Q 

-32m J/4 ,ej/^ ■ p[p P2 ■ p P1 m 2 Q - 32m J/4 ,e J/4 , • p pl p[ ■ p P2 p 3 + 32mj/^£j/^ ■ p P2 p[ ■ p P1 p 3 

-128m J/l> e J/l> • p pi p[ ■ p' 3 p 3 ■ p P2 - Mmj/^ej^, ■ p P2 p[ ■ p 3 p' 3 ■ p pl 

+64m J/4 ,e J/i] ■ p' 3 p[ ■ p P2 p 3 ■ p pl - 64mj /4 ,ej /4 , ■ p\p' 3 ■ p P2 p' 3 ■ p pl 

+32m J / il e J/ ^ ■p' 1 pfp P2 ■ p P1 ] 

+L> M (m 9 , u, v)[96mj/^Sj/^ ■ p P1 p^ 2 m,Q - 32m J/i ,£ J/i , ■ p P2 p p \m 2 Q 
-^mj/^j^pp, ■ p P2 m 2 Q - 32m J/i ,e J/i , ■ p P1 P^ 2 p[ ■ p' 3 - 32m J/i ,e J/i , ■ p P2 p pl p[ ■ p 3 
-32m j 1 j / y, ■ p P1 p'gp\ ■ p P2 + 32m J /^e J/il ■ p 3 p pl p[ ■ p P2 + 32m J/4 ,e J/i , ■ p P2 p£p\ ■ p P1 
+32m J/4 ,e J/i! ■ p' z p p jp\ • P P1 - 32m J/ . l pe J/4 , ■ p Pl p^. 2 p 3 2 + 32m J /^e J/i> -p^p^p' 2 
-96mj/^ej/^ -p pl p^p' 3 -p P2 - 128m J/4 ,£ J/4 , -p P1 p 3 p' 3 ■ p P2 - 32m J/ ^e J/4 , ■ pip^p^ ■ p P2 
-32mj i ^j l ^ ) p\ ■ p pl p 3 ■ p P2 - 32m J/l pe J/l p ■ p P2 pTp 3 ' Ppi - ^ m J/^j/^ ■ P P2 P 3 P 3 • P P1 
-32mj/^£j/^ ■ PiP P2 p 3 ■ p P1 + 64m j/^e J/4 , ■ p' 3 p^ 2 p 3 ■ p P1 + 32m J /^ / ^p[ ■ p P2 p' 3 ■ p pl 
-6Am J/4 ,e^ /4t p 3 ■ p P2 p 3 ■ p pl - 32m J /. i ,e J/i , ■ p' 3 p'iP P2 ■ p P1 + 32m J/i ,e J/4 , ■ p[p 3 p P2 ■ p P1 
+32mj/^ J/ ^p' 1 -p 3 p P2 -p P1 + 32mj/^e^pfp P2 ■ p pl ] 

+D^(m q , u, v)[-6Amj/^£j/^ ■ p P1 p' 3 P P2 + 64mj/^£j/^ ■ p' 3 p^ 2 P pl 

-QQmj/^ej/^ -p pl g^p' 3 ■ p P2 - Umj/^^p^p^- p P2 - 32m J/i ,e J/i , -p^g^p^ ■ p P1 

-32mji^eji^ -p'^p^ ■ p P1 + 64m J/i: e^ /i} p 3 p pl ■ p P2 }} 

If, 1 
D (m q ,u,v) = — Id k 



Cg, U,, V — 

D^rriq.u^) = I d 4 k 

Dn„(m q ,u,v) = t\ / d i k 

ITT 



where p\ = ~P5,P 3 = Pi + P& - P2- 

For amplitudes A lcl and A lc2 , we have 



fc 2 [(fc 


-p 5 ) 2 


-m\]{k 


~P3 


-P5) 2 [(fc- 

kfj. 


hp4 - 


l"P6 


-P2) 2 - 


m o\ 


fc 2 [(fc 




- m 2 ](k 


~P3 


~P 5 ) 2 [(k- 


VPi - 


VPa 


-P2) 2 - 




fc 2 [(fc 


-p 5 ) 2 


- m 2 ](k 


~P3 


-p 5 ) 2 l(k- 


VPi - 


lP§ 


-P2) 2 - 


m Q\ 



A id = C lcl i7 lcl (m gjU;? j)$ JM $p 1 $p 2 
A ic2 = c lc2 H lc2 (m q ,u 1 v)<fj / ^ P1 <fp 2 



£»lcl r^^rparpb^C^r^^rpbrparpC^j 

^jlc2 ___ r£^/rparpb/pc\rj^/rparpbrpc\ 

~ ~ in 2 1 1 

H lcl (m q ,u,v) = -H^(m q ,u,v) = -^d^?^^ 

{E (m q ,u, v)[-32mj/^£j/^ ■ p P1 p 4 ■ P P2 to q ~ 32m J /^£ J /,j 1 ■ p P2 p 4 ■ P P1 m 2 Q 
+96m J/ ^e J/ ^ ■ p' 4 p P2 ■ p P1 m 2 Q - 32m J/ ^£ J/ ^ • p P2 p\ ■ p pl p' 2 ■ p' 4 
-32mj/^ej/^ ■ p pl p' 2 ■ p P2 p[ ■ p' 4 - S&mj/^ej/^ ■ p 4 p[ ■ p pl p' 2 ■ p P2 



-32mj/^ej/^ ■ p P2 p[ ■ p' 4 p 2 
+32m J/ , p Ej / . l p ■p Pl p\ -p' 2 p 4 
+32m J / 1 pSj/. i , -p\p' 2 -p P1 p' 4 
+32mj/ t( ,£j/Tp ■ p 2 p[ ■ p P2 p 4 
-32mj/^£j/^ ■ p' 4 p\ ■ p' 2 p P2 
+32mj/ 1 pSj/ i , ■ p[p 2 ■ p' 4 p P2 



■ p pl - 32m J/4 ,£j/ i , ■ p' 4 p[ ■ p P2 p 2 ■ p pl 

■ p P2 - 32mj/^,£j/^ ■ p' 2 p[ ■ p P1 p' 4 ■ p P2 

■ p P2 + 32m J/i ,£j/ i , ■ p P2 p[ ■ p' 2 p' 4 ■ p pl 

■ p pl - 32m J/4 ,£j /4 , ■ p 4 p 2 ■ p P2 p 4 ■ p P1 

■ p pl + 32m J/i) e J/l p ■ p' 2 p\ ■ p' 4 p P2 ■ p pl 

■ p pl - 32m j i^e j ^ ■ p pl p[ ■ p P2 p 2 ■ p' 4 ] 
+E fl (m g ,u,v)[-32m J /^£j/^ ■p P1 p^ P2 m 1 Q - 32m J / i ,£ J/ ^ ■ p P2 P P1 mQ 
+96m J/i ,e^ /4! p pl ■ p P2 m 2 Q + 32<mj/^£j/^ ■ p P1 p^ 2 p[ ■ p 2 + 32mj/^£j/^ ■ p P2 p%J> x • P2 
-32mjj^eji^ ■ p P1 p^ 2 p[ ■ p 4 - 32mj N ,£ J / 1 p ■ p P2 p^pi ■ p 4 - 32m J/ ^e J/ ^ • p^p^Vi ■ p P2 
-32mj/^£ JN , ■ p P1 p 4 Pi ■ P P2 + Zlmj/^Ej/Tf, ■ p' 2 P P1 p'i ■ P P2 - ^m JN ,£ JN , ■ p' 4 p^ pl p\ ■ p P2 
-32m J/4 ,£j/^ -Pp^Pi -p P1 - 32m J/4 ,£ J/4 , ■ p P2 p 4 p\ -p P1 - 32m J/l> e J/l> • p^Pi ' Ppi 
-mmjj^£ji^ ■ p 4 p^ 2 p[ ■ p P1 - 32m J/l p£ J/i) ■ p P1 p^ 2 p 2 • Pa - 32mj /v ,ej /v , • p^p^p^ ' Pa 
-32m J/4 ,£j/^ ■ p P1 pTP2 ■ Pp2 - 32m J/4 ,£ J/4 , ■ p P1 p 4 p 2 ■ p P2 - 32m J/l> e J/l> • p' x ip% x p' % ■ p P2 
-96mj/^£j/^ ■ p' 4 p^ pl p' 2 ■ P P2 - QGrnj/^jj^Pi ■ p P1 p' 2 ■ p P2 - 32m J/4 ,£ J/4 , ■ p P2 p' 4 p 2 ■ p P1 
-32mj/^£j /4! ■ p P2 p' 4 P 2 ■ P P1 + SZmj/^Ej/^ ■ p'xP^Jp'i ■ p P1 - 32m J/i ,£ J / 4 , ■ p 4 p^ 2 p' 2 ■ p P1 
~32m J/i ,£^ J/i! p' 1 ■ p P2 p 2 ■ p pl + 32m J/i ,£ J/i , ■ p P1 p' 4 p 4 ■ p P2 + 32m J/ . ll) £j /4 , ■ p P1 p' 2 p 4 ■ p P2 
+32m J/4 ,£j/ iJ ■ p'i^p'4 ■ p P2 - 32m J/i ,£j/^ ■ p^p^p^ ■ p P2 - 32mj/^£^p[ ■ p pl p' 4 ■ p P2 
+ 32m J/V £ j/^2 ' PpiP4 ■ Pp2 + 32m J/ ^e J/ ^ • p P2 p' 4 p' 4 ■ p P1 + 32m JN ,£ JN , ■ p P2 p 2 p' 4 ■ p P1 
-32mj/^£j/^ ■ p'iP^ 2 p' 4 ■ p P1 + 32mj/^£j/^ ■ p' 2 p P2 p' 4 ■ p P1 + 32m J /^ J/ ^p' 1 ■ p P2 p' 4 ■ p pl 
-32mj /4 ,£^ J/4t p 2 ■ p P2 p 4 ■ p pl + 32m J/i ,£j/ 1 p ■ p' 2 p' 4 p P2 ■ p P1 - 32mj/^£j/^ ■ p 4 pTp P2 ■ P P1 
+32mj /4 ,£j/^ -p'yP^p^ -p P1 -32m,j/ 1 pEj /l p ■ p 4 p' 2 p P2 -p P1 + 32m J/l> e J/l> -p\p 4 p P2 -p P 
+32m J / 1 p£ J /. i , -p' 2 p 4 p P2 -p P1 -32m J/i ,£^ /%j) p' 2 ■ p[p P2 ■ p pl + 32mj /.^ /%j) p' 4 -p\p P2 ■ p P1 
+32mj /4 ,£^ J/ , ip p' 4 -p' 2 p P2 -p pl ] 

+E llu (m q ,u,v)[~Mm,j N ,£ JN , -p pl p'£p v p2 -&Amjj^£j/^ -p P2 P 4 P P1 
-\2Smj/^£j/^ ■p' 4 p p : 2 p pl - 32mj/^£ J/4 , -p^g^p^ ■ p P2 - 32mj/^£j/^ ■ p P2 g' J ' u p' 1 ■ p P1 
-l2%m J /^ J/%l! p l ' p2 p l -p pl - 32m J / q]; £j/ i) -p^g^p^ ■ p P2 - 128mj/^ /4! p pi p 2 ■ p P2 
^mj/^j/^ ■ p n g^ v p' 2 ■ P P1 + 32TOj /v ,e J/v , • p Pl g^ v p' 4 ■ p P2 + 32m J/ . l p£ J/4 , ■ p P2 g^p' 4 ■ p 
+32m J / i ,£ J/x p -p\g^p P2 ■ p P1 + 32m J/4 ,£ J/i] ■ p' 2 g^ v p P2 ■ p P1 - 32m J/ ^£ J/l p -p' 4 g^p P2 -p 
-Mmj/^e^^p^ -p pl ] 

+E^ e (m q ,u,v)[~32m J/i) £j /4 , ■ p P1 g^p p2 - 32m J/i ,£ J/4 , -p^g^p 9 ^ 
+32m J/ . 4 ,£ t } /4i g" 6 p pl -p P2 - ^rnj/^e^^p 6 ^]} 



E (m q ,u,v) 
■X ( d 4 k 



k 2 [{k+pi) 2 - m 2 Q ][(k+pi -pA-p 5 ) 2 - rn 2 Q ][{k + p 3 ) 2 - m 2 ](k + p x + p 2 - pa ~ Ps) 2 
Ep(m q ,u,v) 

1 f d i k Kl 

in 2 J k 2 [(k +pi) 2 - m 2 Q ][(k+pi - p± - p 5 ) 2 - m 2 Q ][{k + p 3 ) 2 - m 2 ](k + p 1 + p 2 -p± - P5) 2 

' f d 4 k k V k v 

itt 2 J k 2 [(k+pi) 2 - m|][(fc+pi -P4-P5)' 2 - m 2 Q ][{k + p 3 ) 2 - m 2 ]{k + p x + p 2 - Pa - P5) 2 
E^eimq.u.v) = 
^ ^ d^k k^k^kg 



in 2 J k 2 [{k+pi) 2 - m 2 Q ][{k+pi -p 4 -p 5 ) 2 - m 2 Q ][(k+p 3 ) 2 - m 2 ](k + pi + p 2 - p± - Psf 



where p\ = pi 7 p' 2 =px-pi~ P5 7 p' 4 = Pi- 
For amplitudes A 2al and A 2a2 , we have 



A 2ai = c 2al H 2al (m q ,u,v)^ J/ ^ P1 <i>p 2 
A 2a2 = C 2a2 H 2a2 (m qi u,v)<f> J/ ^ P1 <i> P2 
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with 



£j2a2 r^^^/jnarj^brpC^rYY^-parpbrpC^ 

~ ~ in 2 1 1 
H 2al (m qi u,v) = -H 2a2 (m q ,u,v) = — — r 5?(— — ) 3 7 ~-r 

V 9 ^ V q ' (2^) 4ysV 47V c 7 {p 3 + P 5 ) 2 [(pi - P3 - P5) 2 - m 2 Q ] 

{D (m q ,u,v)[96mj N ,e j /ip -p P1 p[ ■ P P2 m 2 Q -32m J/ ^£ J/ ^ ■ p[p P1 ■ p P2 m 2 Q 

-32m J/4 ,£ J/4 , ■ p P2 p[ ■ p P1 m 2 Q - 32m J/i ,£ J/4 , ■ p pl p[ ■ p P2 p 3 + 32m J/4 ,£ J/4 , ■ p P2 p[ ■ p P1 p 3 

-128mj/^ej/^ ■ p P1 p[ ■ p' 3 p' 3 ■ p P2 - 64mj/^sj/^ ■ p P2 p[ ■ p' 3 p 3 ■ p pl 

+64m J/ ^£ J/ ^ ■ p' 3 p[ ■ p P2 p' 3 ■ p pl - 6Amj/^£j/^ ■ p[p 3 ■ p P2 p' 3 ■ p pl 

+32m J/v ,e J/v , • p\p 3 p P2 -p P1 ] 

Dp{m q ,u,v)[-2,2mj/i,Ej/i, ■p P2 ]^ I mg + 96m j/^ej/^ ■ p P1 p P2 m 2 Q 
-32m J/4 ,£^ /4! p pl -p P2 m 2 Q - 32m J/4 ,Ej/^ -p P1 p^ P2 p\ ■ p 3 - 32m J / i ,£j/^ ■ p P2 p^p' x ■ p' 3 

-32m J/v ,£ J/v , • p P1 p' 3 Pi ■ P P2 + 32m.j/,p£j/,p • p'zP p1 p'i • P P2 + 32m J/v ,£j/^ • P P2 P 3 p\ ■ P P1 
+32m J/4 ,e J/4 , ■ p' 3 p^ 2 Pi ■ p P1 - 32mj/^Ej/^ ■ P P1 P^ 2 p'i + 32m J/i ,e J /. ll) ■ p P2 p^ pl p 3 
-96m J/v ,£ J/v , ■p P1 p'{ l p 3 -p P2 - 128mj/^£j/^ ■ p P1 p' 3 p 3 -p P2 - 32m J/i) £ J/l p ■ p\p pl p' 3 -p P2 
-32mj/^e^^i • p P1 p 3 ■ p P2 - 32m J/ ^e J/ ^ ■ P P2 pTp' 3 ' Ppi ~ 64m J/4 ,£ J/4 , ■ p P2 pgp' 3 ■ p P1 
-32m J/4 ,e J/4 , ■ p[p^ 2 p 3 ■ p P1 + QAmj/^Ej/^ ■ p' 3 p^ 2 p 3 ■ p P1 + 32mj /i ,£^ J/ , ifi p' 1 ■ p P2 p 3 ■ p pl 
-64m J/4 ,£^ /4! p 3 ■ p P2 p 3 ■ p pl - 32m J / r4 ,£j/ i , ■ p' 3 pTp P2 • P P1 + 32m J / i ,£j/ i , ■ p[p 3 p P2 ■ p P1 
+32m J/ ^ /i! p' 1 -p' 3 p P2 -p P1 +32m J/i ,£^ J/ , ip p' 3 2 p P2 -p pl ] 

+D fllJ (m q ,u,v)[-64m J/ ^£j / ^ -p P1 p 3 P P2 +64mj/^£j /4 , -p 3 p%j> P2 

-QGmj/^Ej/Tj, ■ p^g^p's, ■ p P2 - 64m j/^^ ■ p" pi p' 3 ■ p P2 - 32m J/4 ,£ J/i! ■ p P2 g^p' 3 ■ p pl 

-32m,ji^£ji^ -p'^p^ -p P1 +64m J /^ /i) p' 3 v p I , 2 -p pl }} 

1 r , 4 , 



D (m q ,u,v) = — Id k 



iir 2 
1 



D,j,(mq,u,v) = -r-j / d A k 



1 


k 2 [(k+ P6 ) 2 


— m 2 ](k - 




+ P6) 2 [(fc- 

kfj. 


I-Pl 


~P3 


~P 5 ) 2 


- m o\ 


k 2 [(k+ P6 ) 2 


— 777,2] (fc _ 


hj74 


+ P6) 2 [(k- 
k^ky 


VPI 


~P3 


~P 5 ) 2 


- m o\ 


k 2 [(k+ P6 ) 2 


— 77*1^] (fc - 


hj74 


fp 6 ) 2 [(fc- 




~P3 


~ P5 ) 2 





D„ v (m q ,u,v) = ^2 / dik T2l71^T^ — ZZwZ I Z rSrTI I Z Z — ZZ2I 



where p[ = p 6 , p' 3 = pi - p 3 - p 5 . 

For amplitudes A 2bl and A 2b2 , we have 



^261 = c 2bl H 2bl (m q ,u,v)<S>, m ,<Z> P1 $ P2 

A 2b2 = c 2b2 H 2b2 (m q ,u,v)^ J/ ^ P1 ^ P2 (29) 
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with 



Qlbl ^j.^^a^brjiC^rY^^rpbrj^arj^C^ 



H (m q ,u,v) = -H ~{m qi u,v) = --—Tg s {——) ■ ■ ^ 2T 

(27r) 4ysV 4iV c y (p 4 + pe) 2 [(p 4 + p 6 - P2) 2 - m 2 Q ) 

{Do{m q ,u,v)[9&mj/^ej/^ • p P2 p[ ■ p P1 m 2 Q - 32m J/l pe J/i) ■ p pi p[ ■ p P2 m 2 Q 
+i2m JN ,£j N , ■ p\p P2 ■ p P1 m,Q + 32m J/v ,e J/v , • p pl p[ ■ p P2 p'£ - 32m j/^ej/^ ■ p P2 p[ ■ p pl p 3 
-eAmj/^ej/^ ■ p pl p[ ■ p' 3 p' 3 ■ p P2 + QAmj/^ej/^, ■ p z p\ ■ p P2 p 3 ■ p P2 
-128m J/ ^e J/v , • p P2 p\ ■ p' 3 p' 3 ■ p pl - Umj/^Ej/^ ■ p[p 3 ■ p P2 p 3 ■ p pl 
+32m J/ ^£j / ^-p' 1 p 3 2 p P2 -p pl ] 

+D^(m qi u,v)[96mj /4 ,ej /4 , ■ P P2 P pl m 2 Q - 32m J/4 ,e.j/^ ■ p P1 p^ 2 m 2 Q 
-Wmj/^jj^p^ ■ p P2 ttiq - 32mj/^£j/^ ■ p P1 p^ 2 p[ ■ p' 3 - 32m J/ ^e J / 4 , ■ p^p^p'x ■ p' 3 
+32m J/ ^e J/ ^ ■ p P1 p 3 p\ ■ p P2 + 32m J/ij £j /i) ■ p' s p^p[ ■ p P2 - 32m J/i ,e J/l p ■ p P2 p' 3 p[ ■ p P1 
+Z2mj/ i ,£j N , ■ p^p^p'x ■ p P1 + 32m J/i ,e J /. ll) ■ P P1 P^ 2 P 3 - 32m J / V ,£j/ V , • P P2 P^ P1 P 3 
-32mj/^£j/^ ■ p P1 p'iP 3 ■ p P2 ~ GAmj/^Sj/^ ■ p P1 p' 3 p 3 ■ p P2 - 32m J/i) ej / . l p ■ p'tf^Pa ■ p P2 
+QAm J / i ,£j/ i , ■ p^p^p^ ■ p P2 + 32m J/4 ,e^ / , ifi p' 1 • p pl p 3 ■ p P2 - 96m J/4 ,ej /4 , ■ p P2 p'iP 3 ■ p P1 
-U&mj/Tpej/^ ■ p P2 p£p' 3 ■ p P1 - 32mj /4 ,Sj /4 , ■ p[p^ 2 p 3 ■ p P1 - 32m J/4 ,e' 1 J/4> p' 1 ■ p P2 p' 3 ■ p pl 
-QAmj/^jj^p'z ■ p P2 p' 3 ■ p pl - 32m J/i) e J/ . l p ■ p' 3 p'iP P2 ■ p P1 + 32m J/i) e J/l p ■ p\p 3 p P2 ■ p P1 
+i2mj / ^ J/ , ip p' l -p' 3 p P2 -p P1 +'&2mj/^ J/ , ip p 3 2 p P2 ■ p pl ] 
+D flu (m q ,u,v)[-64m J /< ll £j/<p ■ p P2 p^ pl p 3 + 64m, 7/ ^e, 7/ ^ ■ p' 3 p P [ 2 p pl 
-Wmj/^Ej/^ ■ p pl g^p' 3 ■ p P2 - 96m J/4 ,e J/4 , ■ p P2 g^p' 3 ■ p P1 - 64m J/4 ,e^ /i: p p2 p 3 ■ p pl 
-32m 7/l/) ej/ l/) • p'^p^ -p P1 + &Am J/i: e^ N> p 3 p pl ■ p P2 }} 

If 1 

D (m q ,u,v) = — I d A k 



in 



D^riq.u^v) = — I d A k 



Dfii/ {rriq , u 1 v) = — I d^k—-— — ^q-j 



k 2 [(k 


-P3) 2 


- m 2 q ] (k 


~P3 


-p 5 ) 2 [(k- 
kp 


hp4 - 


l"P6 


-P2) 2 




k 2 [(k 


-P3) 2 


— m 2 ] (k 


~P3 


-P 5 ) 2 [(k- 
kfj, k v 


hp4 - 




-P2) 2 




k 2 [(k 


~p 3 ) 2 


— m 2 ] (k 


~P3 


-P 5 ) 2 [(k- 


hp 4 - 




-P2) 2 





where p\ = ~P3,P 3 — Pi + P6 - P2- 

For amplitudes A 2cl and A 2c2 , we have 



A 2d = c 2cl i/ 2cl (™ g , U ,z;)$, //V; $p 1 $p 2 

A 2c2 = c 2c2 H 2c2 (m q ,u 1 v)^ JN ^ P1 ^ P2 (31) 
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2 



H™(m q ,u,v) = -H™(m q ,u,v) = j^tii^ 

{E (m q ,u, v)[-Z2mj/^£j/^ ■ p P1 p 4 ■ P P 2 m Q ~ 32mj/^£j/^ ■ p P2 p 4 ■ P P1 m 2 Q 
+96m J/ ^e J/ ^ ■ p' 4 p P2 ■ p P1 m 2 Q - 32m J/l pe J/i) ■ p pl p[ ■ p P2 p 2 ■ p' 4 
-32mj /4 ,ej/^ ■ p P2 p[ ■ p Pl p' 2 • p\ - 32m J/4 ,ej/^ ■ p pl p[ ■ p 4 p' 2 ■ p P2 
-96m J/4 ,ej/^ ■ p' 4 p[ ■ p pl p 2 ■ p P2 - 32m J/4 ,£ J/4 , • p P2 p\ ■ p' A p' 2 ■ p pl 
-mmj/^eji^ ■ p' 4 p' x • p P2 p 2 ■ p pl + 32m J/i) e J/l p ■ p pl p\ ■ p' 2 p' 4 ■ p P2 
+32mj/ i! e J /. 4 , ■ p' 2 p[ ■ p P1 p' 4 ■ p P2 - 32m J/i ,£j/ i , • p[p 2 ■ p P1 p' 4 ■ p P2 
+32m J/4 ,ej/^ ■ p P2 p[ ■ p' 2 p 4 ■ p pl - 32mj/^£j/^ ■ p' 2 p[ ■ p P2 p 4 ■ p P1 
+32m J/i) e J/iJ ■ p\p' 2 ■ p P2 p' 4 ■ p pl - 32m, JN ,£ JN , • p' 4 p\ ■ p' 2 p P2 ■ p pl 
+32m J / 1 pe J / i , ■ p' 2 p[ -p' 4 p P2 ■ p P1 + 32m j/^e jj^, ■ p[p 2 ■ p' A p P2 ■ p pl ] 
+E fl (m g ,u,v)[-32m J /^£j/^ ■p P1 P^ P2 m 2 Q - 32m J/l pE J/ ,p ■ P^P^mg 
+96m J/i ,e^ /4! p pl ■ p P2 m 2 Q + 32<mj/^£j/^ ■ p P1 p^ 2 p[ ■ p 2 + 32mj/^£j/^ ■ p P2 p%J> x ■ p 2 
-32mjj^eji^, ■ p P1 p^ 2 p[ ■ p 4 - 32mj /i ,e J/4 , ■ p^p^p'x ■ p 4 - 32mjj^£jj^ ■ p P1 p 2 p[ ■ p P2 
-Z2mj/^E JN , -p p1 pXp'i -P P2 - 32m, //v ,£ J/v , ■ p' 2 p^ p J>' 4 ■ p P2 - 9&m JN ,£ JN , ■ p' 4 p^ pl p\ ■ p P2 
-32mj/^£j/^ -PpjPsp'x -p P1 -32m J/4 ,e J/4 , ■ p P2 p 4 p\ ■ p P1 + 32m J/l> e J/l> ■ p' 2 p^ 2 p'i ■ P P1 
-2,2mjj^eji^ ■ p 4 p^ 2 p[ ■ p P1 - 32m J/l pe J/i) ■ p P1 p^ 2 p 2 ■ p' 4 - 32m J/i) £ J/lp ■ p P2 ^p' 2 ■ p' 4 
-32m J/4 ,ej/^ ■ p P1 PiP 2 ■ p P2 - 32m J /<pe J /< ll ■ p P1 pfp 2 ■ p P2 + 32m J/l> e J/l> • p' x p% x p' 2 ■ p P2 
-32mj/^ej/^ ■ p' 4 p^ pl p' 2 ■ p P2 - 32m J /^ / ^p 1 ■ p pl p' 2 ■ p P2 - Z2m J/4 ,£ J/4 , ■ p P2 p?p 2 ■ p P1 
-Zlmj/^ej/^ ■ p P2 p' 4 p 2 ■ p P1 - 32m, // ^£ J/ ^ ■ p'^^p^ ■ p P1 - 96™,//^^/^ ■ p 4 p^ 2 p' 2 ■ p P1 
-96mj/^e^pi • p P2 p 2 ■ p pl + 32m J/i ,£ J/i , ■ p P1 p' 4 p 4 ■ p P2 + 32mj/^£j/^ ■ P P1 p' 2 p' A ■ P P2 
-Z2mji^eji^ ■ p' x p% x p\ ■ p P2 + 32m J/4 ,e J/4 , ■ p^p^p^ ■ p P2 + 32m J ^e^p' 1 ■ p P1 p' 4 ■ p P2 
-Z2mj/^ ]/4i p' 2 ■ p pl p A ■ p P2 + 32m J/l pe J/i) ■ p P2 p' 4 p A ■ p P1 + 32m J/ ^£ J/v , • p P2 p 2 p' A ■ p P1 
+32m J/i) £j /i > ■ p'xP^Pi ■ p P1 - 32m J/v ,e J/l> ■ p' 2 p p2 p' 4 ■ p P1 - 32m J/ll ,e^ /it p' 1 ■ p P2 p' 4 ■ p pl 
+32mj /4 ,e^ J/ , 4t p 2 ■ p P2 p 4 ■ p pl + 32m J/v ,e J/ ^ ■ p' 2 p' 4 p P2 ■ p P1 - 32mj/^£j/^ ■ p 4 pTp P2 ■ P P1 
+32m J/4 ,ej/^ -p'lP^p^ ■ p P1 - 32m J/4 ,e J/4 , ■ p 4 p' 2 p P2 -p P1 + 32m J/l> e J/l> -p\p 4 p P2 -p P 
+32m J/i ,£ J /. ll) -p' 2 p 4 p P2 -p P1 - 32m J/v ,£^ /jA P2 ' PiP P2 ' Ppi + 32to J/<A £ j/^P4 'PiPpi ' Ppi 
+32m J/4 ,e^ / , ip p' 4 -p' 2 p P2 ■ p pl ] 

+E llu (m q ,u,v)[~Mm,j N ,E J / xl , -p pl p'£p v p2 -&Amjj^£j/^ -p P2 P 4 P P1 
-USmj/^ej/^, ■p' 4 p p : 2 p pl - 2>2mj/^£ J/4 , -p^g^p'x ■ p P2 - 32mj/^£j/^ -p^g^p'i ■ p P1 
-Yl^mjj^j^p^p'x -p P2 - 32mj/^,£j/^ ■p P1 g' J ' v p 2 ■ p P2 - USmj/^j^p^p^ ■ p pl 
-Wmj/^Ej/^ ■ p P2 g^ v p' 2 ■ p P1 + 32TOj /v ,e J/v , • p Pl g^ v p' 4 ■ p P2 + 32m J/ . l p£ J/4 , ■ p P2 g^p' 4 ■ p 
+32m J / i ,e J / x p -p\g^p P2 -p P1 + 32m J/4 ,e J/i] ■ p^g^p^ ■ p P1 - 32m J/ ^e J/l p -p' 4 g^p P2 -p 
-Mmj/^e^^p^ -p pl ] 

+E^e(m q ,u,v)[~32m J/i) e J/4 , -p^g^p 9 ^ - 32m J/i ,e J/4 , -p^g^p 9 ^ 
+32m J/ . 4 ,£ t } /4i g" 6 p pl -p P2 - ^rnj/.^e^^p 6 ^]} 
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E (m q ,u,v) 
■X ( d 4 k 



k 2 [(k+pi) 2 - m 2 Q ][(k+pi -p 4 -p 5 ) 2 - rri^][(k + p 6 ) 2 - m 2 ](k + pi + p 2 - Pa ~ P5) 2 
Ep(m q ,u,v) 

1 /" d 4 fc K 

in 2 J k 2 [(k +pi) 2 - m 2 Q ][(k+pi -p 4 -p 5 ) 2 - "^][(fc + P6) 2 - "^](Hpi + p 2 - Pi - P5) 2 

E ftl/ (m q ,v,,v) 

1 / d 4 fc fc M fc ^ 

i7r 2 7 fc 2 [(fc +pi) 2 - m^][(ft +p x -P4-P5) 2 - m 2 Q ][{k + p 6 ) 2 - m 2 ]{k + p x + p 2 - Pa - p 5 ) 2 

^ f d^k kfj,k v kff 
in 2 J k 2 [(k +P1) 2 - m 2 Q ][(k + p 1 - p 4 - p 5 ) 2 - m 2 Q ] [{k + p 6 ) 2 - m 2 ](k + pi + p 2 - Pa - Ps) 2 

where p[ = pi,p 2 = Pi - Pa - P5,P 4 = Pe- 



(33) 



2. The case of J/V> - 
For amplitude A lala , 



VP 

we have 



A 1 



C 



lala Tjlala 



{m q ,u,v)<$>jj^ v <&i 



(34) 



with 



H lala (m q ,U,v) 



Tr (T" 2 ? 1 T' c )Tr (T" 3 "?^ T" c ) 

(2ir)* gs[ 4Nc ' ( P3 + P5 ) 2 [(pi -ps - P5 ) 2 - m 2 ] 
m Q m, z {L>o(m g ,u,w)e Aira/3 [-64e(; /i> £*,»^3 -p^ - 64e^ /v ,e; a ^ /v ,^p3 • p p 
+64 £ ^<«p^p p p 3 -p v + 64 £ ^<«p 3 "p p p J/v , -p v - 64 £ ^<%^3V ' Pp] 

+D v (m q ,u,v)e^ a p(32e$ / ^e™pPp J/ i, -p p - 32e^ /i) e* v a p l3 p pj / ^-p v 
-32e^e*/p<}^p v -p P )}} 



(35) 



where p[ =pa,p' 3 = Pi - 
For amplitude ^4 lalf> , 



"P3 ~P5- 

we have 



.4 



lalb 



C lal 'ff lall (m S) tt,B)$ w *v*p 



(36) 



with 



C 



lalft 



Tr(T ia T'^'J lc )Xr('J ia T^2^ c ) 



ff lolt, (m 9 ,ii,u) 



(27T)4 * s * 4N C ' (p 3 + p 5 ) 2 [(p x - p 3 - p 5 ) 2 _ m 2 ] 

m Q m 9 {Oo(m„ u, u ^^[^^^p^Pa ■ p J/v , - 64 £ ^ £ ; t y; / ^p3 • p p 

+64^^ e ;«p^j/ 3 -p v - 64^4^^ -p p + 64e(;^ £ ;' 3 p« /v ,P3^ v -p, 



[D e ( mq ,u,v)e^(96e^py v p e J/l , - S&^O^p" - 32^^/^) 



a e 



+D v {m q ,u, v)e^ a f 3 {32e f * J/ij e* v a pPpj / ^ ■ p p - 32^8 -p^pj/^ ■ p v 
+32e^e*J>p a J/4 P v -p P )}} 



(37) 



where p[ = pa.p'z =Pi~P3- Ps- 



For amplitude A la2a , we have 

A ia2a = c la2a H la2a (m q ,u,v)<f> J/4 ,<I>v<S>p 

with 

ffla2a/ „, „\ _ r fi ( 

h (m g ,u,v) - 77mSs(.-n^~) 



(2tt) 4 ^ V AN C > (p 3 + P5 )2 [(p x - p 3 - p 5 )2 - m 2 ] 

m Q m,{D (m„ u, w)e (tit/Q ^[64e^ /v ,£*^^p3 • p J/4 , + 64e^< a ]» j/^Ps • p p 

+D"(m qi u,v)e^ al3 (-32e^ /tlj e* a pPpj/^ ■ p p - 32e^e*> P pj/^ • p v 
-32 £ ^4%^ v -p P )]} 



where p[ = p 4 ,p' 3 = Px~P3~ P5- 
For amplitude A la2b , we have 



A ia2b = c la2b H la2b {m q ,u,v)^j,^ v ^p 

with 

H™\m q ,u,v) = -gl^(-i-)3 1 



(2tt) 4 y * v 4iV c ; (p 3 + P5 ) 2 [(p! - p 3 - P 5 ) 2 - m|] 
m Q m q {Do(m q ,u, v)e fil/a ^[6Ae^e* v l 'p^p^p' 3 ■ p J/4 , - 64e^e* a ]j^j^ ■ p p 
-64^^4>^p^3 • p v - Me^^ptp^pj,^ • P P + Me^efp^tfpv ■ p P ] 
+ [D e {m q , u, v)e^(-64e^s* v v py^ ~ Z2e» ^p^ ^ v p p + Zle^^p^^p 
+D v (m q , u, v)e^ afi (32e^ /i) £* v a p^p J/i) ■ p p + 32£^ /v ,£* a p P p J/ , / , • p v 



where p[ = p 4 , p' 3 = Pi - P3 - P5 ■ 
For amplitude A lbla , we have 



A ibia = c lbla H lbla {m q ,u,v)^j/^ v ^p 

with 

^jlbla r^^^rj^arjibrpC^rj^^rjibrparpC^ 

iTT 2 a, I 



H (m q ,u,v) - -j——-7g s \~TTi~) 



(2vy ysK ANc> (P4+P6) 2 [(P4+ P6 -P2) 2 - m 2 Q ] 
m Q m q {D {m q ,u,v)e flvaP [64:e' 1 J/i> el a pj /i) p^p / 3 -p p + 64,e^ j/ip e* v a pj /ip p p p 3 ■ p v 
+MEV ]/ ^: a ptp P Pj H ,-p v -^ J/ ^p'»p a ] ^p v -p P ] 

+[Dg(m q , u, W )e^,(32e^<>>e^ - ^j/^TPj/M + ^K/^Pj/M 
+D' y (m q ,u,v)s fll , al 3(-32e^ /i! el a pPpj /4 , ■ p p - 32e^e* Q p^p J/v , ■ p v 
-S2e^efp a J/i!Pv -p P )}} 



where p[ = -p5,p' 3 = Pi + P6 - P2- 
For amplitude A lblb , we have 



A ibib = c lblb H lblb {m q ,u,v)^j/^ v ^ P 
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with 

H^\m q ,u,v) = -gjflJC 1 )' 



(2tt)4 ys * AN C > ( P 4+P6) 2 [(P4+P6-P2) 2 - m%] 
m Q m q {Do(m q ,u,v)e^ a i3[-64:S^ /i! e* v a pj / ^,p' 3 ■ p p - Q^^e^pj^p^ ■ p v 
-64e^e* v a p'fpPp J/4 , ■ Pp + 64e^e*J 3 p^p a J/ ^p v - Pp ] 

+D u (m q ,u,v)e^ a f 3 (32e^ /ip sl a p l lpj /4 , ■ p p + 32e^,e* a p£p, //v , ■ p v 

+WZj/^Pj/TpPv-Pr)}} (45) 



where p[ = -p 5 ,p' 3 = Pi + p 6 - p 2 - 
For amplitude A lb2a , we have 



A uua = C' bZa H' bZa {m qi Uiv)<$>jj^ v <S>p (46) 

with 

^j\b2a r£y^rpap^bpc^rY^^rpapbpc^ 

jjU2a( \ 6/ 1 \3 1 

[m q ,u,v) - {27r)i 9s[ 4Nc ) {p4+p6)2[(p4+p6 _ p2) 2_ m 2 Q] 

m Q m q {Do(m q ,u,v)e^ a p[-64e^e* v a p v J/ ^p 3 ■ p p + 64£^£* a p^ /v ,p p p^ .p v 

+[I> fl (m g ,« > W )£^ Q/3 (-32£^£>>£p^ + S^O^^t + 32e W<>J//* 
+D v (m q , u, v)e^ afi (32e^ /i) £* v a p^p J/i) ■ p p - 32^ /v ,£* a p^p J/v , • p v 

^K/^P^PvpM (47) 

where p[ = -p5,p' 3 = P4 + P6 - P2- 
For amplitude „4 1&2b , we have 



^iwt = c^H^tm^u^v^j/^vSp (48) 



with 

■yl626 ITU /rriarribr 



ulb2b( \ ^ 6/ \ \3 1 

1 9 ' ' J " (2^^4iV c j ( P4+P6 )2[ (p4+p6 _p 2 )2_ ro 2 



m Q m 9 {Do(m 9 ,u,w)e^ C[) 3[-64£^£; Q p^ / ^p; i • p p + 64e^ J/ ^E* v a p"j /i ,p p p 3 ■ p v 
-64e^e* v a p 3 ^ vPj/4 ,-p p + 64e^£*/p^ /v ,p v -p p ] 

+[D e K,u ) «)£ M ^(324 /i/ ,£»^pS / ^ - 32^e^/^4 - 32£^£>^X/J 
+D v (m q ,u,v)E^ a(3 (32E^ /tlj s* v a pPpj /4 , ■ p p - 32£'j / , i/ ,£*>^p, // ^ -p v 



+32£^£*/p^ v ,. Pp )]} (49) 



where pi = -p 5 ,p' 3 = p 4 + Pe - P2- 
For amplitude A lcla , we have 



41cm = c , ^^ij^^( TOgiU , w )$ J/v ,$ y $ p (50) 



with 



m Q m,{£? (m„ii,i;)e MVQ/ 3[-32e^e*,''p°p^p^ ■ - 32e^e; a p^ /v ^^ -p p - 32e f } / ^e* v a p u JN ,p l3 p p' 2 ■ p 
-32s^s^p a M ■ pj/^ - 32 S » /ilS * v a p»j / y v p' 1 ■ Pp - 32e^Oj/^Pi ' Pv 
+32e^e* v a pT4pj^ ■ Pp + 32e^e* v a p?pPpj^ ■ p P - 32^<V 2 ^Pj /V , ■ P v 

-^"j/^TPlPpPJ/i* ■ Pv ~ ^J/^PiPj/^Pp 'Pv - ^J/^PlPj/^Pp 'Pv] 

+[£i*K,K,<;) e ^(-64 £ ^ 

+EZ(m q ,u,v)£^ a p(-64£^ £ * v a pPp J/ ^.p p + 64£^ -p p )}} 



where p[ =Pi,p' 2 = Pi -P4~P5,P4 = Pz- 
For amplitude A lclb , we have 



with 



A lclb = C lcU H lcl \m q ,u,v)<5>j N ,<$> v $p 



fl^m,,^) = -i^(^)3 1 

(27r) 4 4iV c (P4+P5) 

m Q m q {E (m q , u, v)£^ a0 [-32£ f } /ip £* v p^ rP ' 2 ■ p J/4> + 32e^ / ^Pj/^p^ ■ P P + 32s^ /ii s* v a p v J/i! p p p 2 ■ p 

-32 £ ^ £ »vPi ' Pj/i> + ZK/^*vPj/y v Pi ■ Pp + 32e ' Pv 

+32e^e* v y' 2 ^ v p J/4 , ■ Pp + 32e^e* v a p'fpPpw ' Pp ~ ^j/^Tp'ffejf* ' Pv 
-32e v J/ ^p l ?pipj^-p v +32e»j / ^p'£p a J/ ^p P ■ p v + 32e» J/i ,e*J } p?p'}^p p ■ p v ] 
^E ie {m q ,u,v)e^ a p{-&^ J/ ^ P yy m 

+EZ(m q ,u,v)£ flua p(-64£V /ij£ * v a pPpj /i! . Pp +64£^ /ii £* v a p p p J/i) ■ p v - Ue^^efp^^Py -p p )}} 



where p[ = p\,p 2 =P\-p±- P5,p' 4 = P3- 
For amplitude A lc2a , we have 



with 



A vcia = c lcZa H lcZa (m q ,u,v)^j /4 ,^ v ^p 

m Q m q {E (m q ,u, u)e^ Q/3 [32e^£^p^ ■ p J/v , + 32e^ /j/; e^ /v ^^ ■ p p - 32e^ / , e;>^ / , iA p^^ • p v 

+3K/^7pM ■ PJ/O + 32e V> ' Pp ~ ^j/^TPj/^Pi ' Pv 

-M£j/^7pTp v PJ/tP -Pp - ^j/^TpTp v Pj/^-P P - ^j/^TpfppPj/^ " Pv 
-^j/^TpTPpPj/^I> 'Pv - 32ej / ^e;V 2 AI Pj/^ 'Pv ~ ^j/^PiPj/^Pp ' Pv] 
+ [Eie{m q ,u,v)s^ a p(64:Sj^e* v l/ p pP P P j/ 1 p + ^j/^^Pj/^PyPp ~ ^j/^TPj/ipPpPt^ 
+EZ(m q , u, v)e^ a p(64e<; /ip e* v a pPpj /4 , ■ Pp + 64e^e* a p£p J/v , ■ p v + &Ae^efp a j^p v ■ p P )}} 
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where p[ = p 1 ,p' 2 =p 1 -p i - p 5l p' 4 = p 3 . 
For amplitude A lc2b , we have 



A ic2b = C lc2b H lc2b (m q ,u,v)<5>j N ,<$> v §p (56) 



with 



£iic2b Xr(T ,a T'^T' c )Xr ^2 na T'^T ,c ) 

H^ b (m g ,U,v) = J^gH-^f 1 
(27r) 4 4iV c 

mQmg{-B (m g) u,u)ff^ Q:/ 3[32£^ /i/) e v p p p' J v p 2 ■ pj/^, + 32e l } /ip e v Pj^P P v P 2 -P P - ^j/^ v Pj/^P r P2 • Pv 
+32^ /V)£ »^' 1 • p J/4 , + 326^8^^ ■ Pp - • Pv 

+32^0^^ -p P - ^j^p^yPpj^-p.+^^p'^pj^-p, 
+2te'j H tf t p'£if p pj, il ■p v +32E^s* v %^ /4 p P -p v +32s^el%^p P -p v ] 
^E ie (m q ,u,v)e^(64s^p a p p^ 

+E%(m q ,u,v)£ fJ , ua p(-64:e'$ /ii e* v a pPpj / Tp ■ p p - 64e^£* Q ^p, 7/ ^ • p v - QAe tl J/i) e* v p p a J/lp p v -p p )]} 

where p[ = pi,p 2 =Px~P4 -P5,Pi = P3- 
For amplitude A 2ala , we have 



(57) 



A iava = C Zala H Zala (m q ,u,v)$j]^v$p (58) 

with 

1 „ 1 



H 2ala (m q ,u,v) = 7 __ r5 6 (__)3 

V q , ) ( 27r )4 y ^4iV c ; (p 3 + P 5 ) 2 [(pi - P3 - P5) 2 - m 2 Q ] 

mQm q {Do(rn q ,u,v)e^ a ( 3 [-6A£ f } /ip £* v a p u j /i ,p l l_p' 3 ■ p p + GAe^^p^^p^ ■ p v 

+^jl^X a p'£p l3 p P m ,-p v ~ Me^efpZp^Pr -P v \ 
^D e {m q ,u,v)e^ aP {-Z2e^p^ J/t 

+D"(m q ,u,v)e^ a p(32e p J/ip el a p1p J/il ■ p p - 32e^ /i/ e* Q p^ J/ ^ ■ p v 

-i2e^^p a j^p v -p P )]} (59) 

where p[ = p 6 ,p 3 = pi - p 3 - p 5 . 
For amplitude A 2alb , we have 



^2aib = c 2alb H 2alb (m q ,u ,v)^j/^ v ^ P (60) 

with 



( 27r )4^V 4 iV C y (p 3 +p 5 ) 2 [( Pl -p 3 -P5) 2 - 

m«m,{i)o(ffl g , u, v)e^ a p[~64e^ J/4! e* a p u J/i ,p^p' 3 ■ Pp + Q&e> 1 J/ ^e* v a p u j /%ll pP p p 3 ■ p v 

+ [D e (m q ,u,v)e^ aP (32e^ P yy j/i; - 32 £ ^ £ -j^/ p " ^j/^TPj/M) 
+D"(m g ,u, v)e tlva0 (32e^e* v a p^p J/ ^ ■ Pp - 32e^e* v a p p p J/i , ■ p v 

+32e^<%^- Pp )]} (61) 



where p[ = p 6 ,p' 3 = pi-p 3 - P5- 
For amplitude A 2a2a , we have 

A 2a2a = C 2a2a H 2a2a (m q ,U,v)$j /i ,$ V $ P 

with 

fj2a2a/ \ * 71 " 2 6/ ^ \3 ^ 

if (m q ,u,v) = — — rrffs ( TTT - ) 



(27T)4 ysV 4iV C ; (p 3 +p 5 ) 2 [( Pl -P3-P5) 2 -m|] 

m Q m q {D (m q , u, v)e llua i 3 [Me^ / ^e* v a p v ]/i) p fi v p' z ■ p p + 64e^e* a p^p£p£ ■ Py 
+[iMm,,« > «) e |. ra/J (32 e V^ 

+L> l '(m (? , u, u)e M ^(-32£^ /v ,£* a p^p^ • p p - 32£^ /v ,£*>^ J/ , / , • p v 
-32e^p a J/i)Pv .p P )}} 



where pi = p 6 , p 3 = pi - p 3 - p 5 . 
For amplitude A 2a2b , we have 



A 2a2b = C 2a2b H 2a2b {m qi u,v)^j/^ v ^ P 

with 

rrlalbi \ _ _ n 6 ( _}_\3 \ 

{ " ' } ~ (2tt)^W (p 3 +p 5 ) 2 [(p 1 -p 3 -p 5 ) 2 -m 2 Q } 

m Q m q {Do(m q ,u,v)e^ a p[-64e^e* v a p v j / ^p 3 ■ p p - 64e>* J/ip e* a p' / j /4 ,p P 'p 3 ■ p v 
-eie^e^pfpPpj^ ■ p v + Me^e*M v p a j^P P ■ P v ] 

+ [D 6 (m q ,u, v)e^ a0 (-32s^ P yy j/4 , - ^^fy^Pfi + 32£^£>^pX 
+D v (m q , u, v)e flvaj3 (32e^£* v a p^p J/i , ■ p p + 32e^ /4! e* a p p p p J/ ^ ■ p v 
+S2e^e*/pJ^p v -p p )}} 

where p[ = p 6 ,p 3 = pi-p 3 - P5- 
For amplitude A 2bla , we have 

A 2bla = C 2bla H 2bla {m q ,U,v)$j M $ v $ P 

with 

^j2bla ^Yy^pd'J^brpC^rY^^rparpbrjiC'^ 

H [m q ,u,v) - - g s ( 



(2tt) 4 ys y AN C ' (P4+P 6 ) 2 [(P4+P 6 -P2) 2 - m 2 Q ] 
m Q m q {D (m q , u, v)e t , vaP [64 : e^ /i> e* v l/ p p p^p 3 ■ p J/4 , + ME^^pj^p^p's ■ p P 
+Me^ / y v a p" J/ y p p' 3 ■ Pv + ME^e^pPpj^ ■ p v - QAe^ J/ ^fp a J/ ^p v ■ p P ] 
+[D e (m q , u, v)e ftva0 (-32e^ P y v p e JM - 32^-f^p* + 32^0^^ 
+D v (m q ,u,v)8^ a t,(-328^ /ip sl a pPp J/i! -p p - 32e^ /i> £* Q ^p, 7/v , ■ p v 
-32sys*/py Pv -p P )]} 



where p[ = -p 3 ,p' 3 =Pi+P&-p2- 
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For amplitude A 2blb , we have 

A 2blb = C 2blb H 2blb {m qi u,v)<S>,j/^ v ®p (68) 



with 



H (m q ,u,v) - - , 4 9 s {-n^~) 



(2tt)4 ys v 4N C ' (P4+P6) 2 [(P4+ P6 -P2) 2 - m 2 Q ] 
m Q m q {D (m q , u, v)e tlva p[QAe' 1 J/%ji e^p a p p v p' 3 ■ p J/4 , - 64e^e*°p^,p£p3 • p p 
-64e^e* v a p^p p p' 3 -p v - QAs^e^pfp^pj/^ ■ p v + 64 £ ^<^f#p v -p p ] 
+[D fl (m„ u, «)^(-9^0«p^p^ - 32e^4>^p P + 32^4>^pX ) 
+D h '(m q ,u,v)e ^ a p{32e^ /t/j £* v a p v Pj /V '-Pp + 32e j/^r?p^/</' 'Pv 

+32e^<.% v ^.p p )]} (69) 

where pi = -p 3 ,P3 = P4 + Pe - Pi- 
For amplitude „4 2h2a , we have 

A 2b2a = C 2b2a H 2b2a (m qi u,v)<5>j/^ v § P (70) 

with 

^j2b2a ^^^jiarpbrpc^rY^^rpbrparpc^ 

Trivial m ,, „\ _ _^_^6/_!_\3 1 
# [m q ,U,V) - TTTTiS'^TT^) 



(2 7 r) 4 '^ v 47V c / (p 4 + p 6 ) 2 [(p 4 + Pe -P2) 2 - m 2 Q ] 
mQm q {D (m q ,u,v)e^ a ^[-64e^e* v l 'p^p^p 3 ■ p J/4 , - QAs^^p^^plp^ ■ p p 
+64 £ ^4>^p P p 3 ■ Pv + ^j/^p'^pj/^ ■ p v - Me^efp^p'^ ■ Pp ] 
+ [D e (m qi u, v)e, va ^ ]/ ^pyP v p e J/ ^ + We^p^pi + 32e^<>^p p pt) 
+D u (m q ,u,v)e fn ; a {3{32e'" J/ip e* v a pPpj /i , ■ p p - 32e^e* a j^pj^ • p v 

-ZK/^fPj/^VPr)}} (71) 



where p[ = - P 3,P 3 = Pi + Pe - Pi- 
For amplitude A 2b2b , we have 



A ibib = c MM, H M *>(m q ,u,v)$ J/i! $ v $p (72) 

with 

£l2b2b ^^^parpbrpC\rT\^^rj^brparjiC^ 



H 2b2b (m u v) - l " a 6 ( 1 ) 3 

{ m q ,u,v) - (27r) ,9 s [ mc ) ( p 4+ p 6 )2 [(p4+p6 _ p2) 2_ ra 2 



m Q m q {D (m q , u, v)e l _ iua! 3[~64e^ ^e" 'p p p^p 3 ■ pj/^p - QAe^^e^p^j^p^ ■ p p 
+64e^e* v a p^ p p' 3 ■ p v - 64^^<>^ p p J/v , • Pv + die^efp^p'^ ■ p P ] 

+[De(m gi u,v)e MVafi {96^ J/ ^p^y j ^ - 32 e ^ e *«#/X*£ " ^j/^Pj/M) 
+D'- , (m q ,u,v)e^ a p(32e^ /i> e* v a pPp J/4 , ■ p p - 32£^e* Q p p p J/ , / , ■ p v 

+32 £ ^ £ fp^p v ,.p p )]} (73) 



where p[ = -p3,p' 3 = p 4 + p 6 - Pi- 
For amplitude A 2cla , we have 



A iaa = c 2cia H 2cia {m qi u,v)^j/^ v ^ P (74) 



with 



m Q m,j{i;o(™ 9 ,u,«)£^/3[- 32£ j/^ £ t»v-P2 ■ Pj/v ~ 32e j/</>< a Pj/V<PvP2 'Pp - ^j/^TP'j/i-PpP^P 

-^"j/^TPiPpPJ/i* ■ Pv ~ ^ eV j/^ £ *vP2P a j/i,Pp -Pv ~ Mej/^PiPj/^Pp -Pv] 
+[£i*K,K,<;) e ^(-64 £ ^ 

+EZ(m q ,u,v)£^ a p(-64£^ £ * v a pPp J/ ^.p p + 64£^ -p p )}} 



where p[ =Pi,p' 2 = Pi ~P4~P5,P4 = Pe- 
For amplitude A 2clb , we have 



with 



A 2clb = C 2clb H 2clb {m q ,u,v)<5>j N ,<S> v $ P 



m Q m q {E (m q ,u, v)e^ a0 [-32s^ ^e*^ ' PpP ^p 2 ■ Pj/^ + 32e^ / ^e* v a p v JN ,p f3 v p l 2 ■ p p + 32s* /ii s* v a p v J/i! p P 'p 2 ■ p 

-32e ' + 32£ ' Pp + ^j/^vP V J,y P P\ ■ Pv 

+326^0^^ • Pp + 32e^el a pTp v PJ/4> ■ Pp ~ ^j/^T^Pj/i, ■ Pv 
-32e v J/ ^p l rpipj^-p v +32e v j / ^p^p a J/ ^p P ■ p v + 32e» J/i ,e*J } p?p'}^p p ■ p v ] 
^E ie {m q ,u,v)e^ aP {-&^ J/ ^pyy m 

+E^(m q ,u,v)£ flua ^(-64£^ /ij £* v a p^pj /i , -p p +64£^ /i) £* v a p l3 p pj /i) ■ p v - Ue^^efp^^Py -p p )}} 



where p[ =Pi,p' 2 = pi -P4~P5,P4 = Pe- 
For amplitude A 2c2a , we have 



with 



A *ua = C 2cZa H Zc2a {m q , Ul v)<$>j /4 ,<$> v <$> P 



H 2c2a (m q ,u,v) = --^- g e s (—) 3 - 

V «, , ) ( 27r )4 y ^47V c y (P4+P5) 2 

m Q m q {E o (m q ,u,v)e^ a0 [32e , } /ip e*^ P y^p 2 ■ Pj/i , + 32s» /i) e* v a p v J/il pPp 2 ■ p p - 32s» /^Ppp'i ■ P v 

+32^/^»^i ' PJ/O + S^Yj/ZvPl • Pp - We^p^pPp', ■ Pv 

-M£j/^7pTp v PJ/tP -Pp - ^j/^TpTp v Pj/^-P P - ^j/^TpfppPj/^ " Pv 

-^j/^TPiPpPJ/^ 'Pv - ^j/^pTpj/^Pp 'Pv - Ztej/^PiPj/ipPp 'Pv] 

+ [Eie{m q ,u,v)s^ a p(64:Sj^e* v l/ p p : pPpj/ 1 p + ^j/^^Pj/^PyPp ~ ^j/^^Pj/^PpPv) 

+EZ(m q , u, v)e^ a p(64e<; /ip e* v a pPpj /4 , ■ p P + 64e^e* a p£p J/v , ■ p v + &Ae^efp a j^p v ■ p P )}} 
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where p[ = p 1 ,p' 2 =p 1 -p i - p 5l p' 4 = p e . 
For amplitude A 2c2b , we have 



with 



A 2c2b = C 2c2b H 2c2b {m q ,u,v)^,j/^ v ^p (80) 



H 2c2b (m q , u, v) = --^gH-^-f 1 

(27r) 4 4N c [Pi + PbY 

mQrn q {E o (m q ,u,v)s tlva0 [32e^ /jjj s* v u p^p' 2 ■ Pj/i , + 32e^ /ip s* v a p" JN ,pPp' 2 ■ p p - 32s* l ^e*™p v j j ^pP p p 2 ■ p v 

+32£ ' *W + 32e J/^vPj/^tPi ■ Pp - ^j/^TPj/^A ■ Pv 

+W£j,^7pTp v PJ/O ■ P, + ^j/^Tp'Mpj/^ ■ Pp + 32e VO^p^M ' Pv 
+32e^e* v a pTp P Pj/^ ■p v +32e" J/ ^fp^p a J ^p P -p v +32e v J/ ^p l i l p a m p P -p v ] 
+[E ie (m q , u, v)s^ a/3 (64e^ P y v p e J/il + 6^s* v a p^ P y p - ^^p^p^) 
+E'((m q ,u,v)e^ a0 {-Me^ J/ ^e* v a p p v pj/^ ■ p p - QAe^^pPpj/^ ■ p v - 64e^ J/ip e*fp a J/4 ,p v ■ p p )}} 



(81) 



where p[ = px , p' 2 = p x - p 4 - p 5 , Pa 
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